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Abstract. We consider an exactly solvable discretisation of the radial Schrodinger 
equation of the hydrogen atom with 1 = 0. We first examine direct solutions 
of the finite difference equation and remark that the solutions can be analytically 
continued entire functions. A recursive expression for the coefficients in the solution 
is obtained. The next step is to identify the related three-term recursion relation for 
Pollaczek polynomials. One-to-one correspondence between the spectral and position 
representations facilitates the evaluation of Pollaczek polynomials corresponding to the 
discrete spectrum. Finally, we obtain two alternative and explicit expressions for the 
solutions of the original difference equation. 
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1. Introduction 

The general solution of the Schrodinger equation, the corresponding energy levels 
and resulting atomic shell model are taught to all students of quantum mechanics. 
Experimental and theoretical studies of Rydberg atoms with a single electron have 
illuminated the wave functions and the classical limit of quantum mechanics. 

Here, we will concentrate on a specific aspect of theoretical studies, namely 
discretisations and orthogonal polynomials. Orthogonal polynomials, namely Laguerre 
polynomials, are already present in the solution of the Schrodinger equation. The subject 
of orthogonal polynomials related to the hydrogen atom was reviewed in 1991 by Dehesa 
et al. in Ref. pp. Apart from polynomials present in the solutions, different kinds of 
discretisations induce further orthogonal polynomials. The most obvious ones arise from 
direct discretisations of the Schrodinger equation. Unfortunately, the explicit form of 
these polynomials is not usually known. The continuum states and L 2 discretisations of 
the continuum have been studied already in the 1970's [21 IS]- Recently, a corresponding 
solution for the Dirac-Coulomb problem was presented in Ref. jlj. In the context of 
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condensed matter physics discretised Schrddinger equations can also be used as simple 
models of nanoscale systems or related to the tight-binding approximation j^j. 

In this paper, we show that the I = states of the hydrogen atom can be 
exactly and explicitly obtained for the symmetric discretisation of the second-order 
derivative. First this is done by inserting an explicit ansatz in to the difference equation. 
Later on, we relate the difference equation to the three-term recursion relation for 
Pollaczek polynomials with specific parameters. Finally, a surprisingly simple and 
explicit expression is obtained for Pollaczek polynomials and the corresponding solutions 
of the discretised Schrddinger equation. 



2. Initial steps 

Let us consider the radial Schrddinger equation for the hydrogen atom for the hydrogen 
atom, i.e. 

_ h 2 R"(r) _ h 2 R'(r) _ e 2 R(r) + h 2 l(l + l)R(r) = 
2m mr Ane^r 2mr 2 

In so-called natural units and for I = the equation simplifies to 
u"(r) u(r) , . 

- ~y - -y = ( 2 ) 

where u(r) = rR(r). The simple eigenvalues and well-known solutions are expressible 
in terms of associated Laguerre polynomials of the first kind 

«n(r) = rLi(2r/^^^^ E n = (3) 

where n = 1,2,.... Next we discretise this equation using the symmetric second-order 
difference and obtain a finite-difference equation 

- u(r-t)-Mr)+u(r + 5) _ «(r) = (4) 

2o 2 r 

In the following section, we lift the restrictions that r and 5 must lie on the positive 
axis and allow complex values for both. Of course, any solution of Eq. (jlj) multiplied by 
an arbitrary function with period 5 is still a solution. Nevertheless, we will concentrate 
on solutions and eigenvalues that tend to corresponding classical solutions in the limit 
5^0. 



3. Solutions in coordinate representation 

For nonzero values of 5 we note that functions u(r) are also solutions of 

u(r-5)/2 + u(r + 5)/2 + 5 2 u(r)/r = fiu(r), (5) 

where \x = —5 2 E + 1. This problem has been studied by Berezin in the case of purely 
imaginary 5 in Ref. 0. We discovered the existence of an explicit solution in Ref. [Hj. 
Here a more transparent and instructive derivation is given and we are able to obtain a 
new, recursive formula for arbitrary terms in the solution. 
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Let us insert an ansatz 

M (r) = e^f> fc r fc , (6) 
fc=i 

into Eq. (jSJ) and assume 7^ 0, which obviously corresponds to the solution u n (r) in 
Eq. (jSJ). Now the equations must hold identically in r so each coefficient of r J must 
vanish. This yields 

E 7i, fc « fc = 0, 7i+i,fc := <^~ j ( * ) (e* + {-!)*-* e^)/* - ,ryw + * 2 7*j + i (7) 

fc=max(l,j — 1) V 3 f 

for j = 1, . . . , n + 1. The (n + l)th equation requires that a n [cosh(<5/3) — /i] = 0, so we 
find /i = cosh(5/3). Next equation is then simplified to a n [n5 _1 sinh(<5/3) + 1] = 0. These 
equations now yield 

= v/l + (5/n) 2 , = -arsinh(5/n)/5 (8) 

in agreement with Ref. jHj- The classical limits E n (8) — > — l/2n 2 and /3 n (5) — > — 1/n are 
also satisfied. The remaining n—1 equations can be used in order to solve the constants 
«fc with k = 1, . . . , n — 1. Note that the present derivation is both simpler and more 
exhausting than the previous one, which was based on an intelligent guess concerning 
the identity of terms in the power series used. All eigenvalues and constants in the 
exponential part are identical to those given in Eq. (jSJ). 
The general solution to Eq. (J3j) now becomes 

= (E ^ «i n ' 5) r fc ) exp(-r arsinh(5/n)/5), := - (jlj)- ( 9 ) 

The coefficients {a^'' 5 ''} are of the form 

«i n i = (i + *7- 2 ) (fe - 2fc ' )/2 E 41> 2m , (io) 



where atjj. q = 1 and fc' = |_ A;/ 2 J , i.e. fc/2 if A; is even and (k — 1)/2 if fc is odd. In addition, 
we find 



(n) -2m 

a n-k,m ~ n 



[k/2\\ P(2m-l)(n-k)\ 
m J (n-k + 2m- 1)! ' 

where P(2m — 1) is a polynomial of order 2m — 1 such that the coefficient of x 2ra ~ x is 
equal to unity. The general form of the leading terms is given by 

«i n ' 5) = 1 (12) 
= v/l + 57^ (13) 
(3n - 1)(5 2 
3n 2 (n — 1) 



(n,s) 1 , (3n-l)£ 2 



a <»4».^A?( 1 + -^) (15) 

(n,g) 2(n - 1)£ 2 (15n 3 - 30n 2 + 5n + 2)£ 4 

an - 4 n 2 (n-3) 15n 4 (n-l)(n-2)(n-3) 1 J 
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Here n is an arbitrary state index, so once ctn-1 nas been obtained, we have exact 
expressions for k + 1 leading polynomial terms in any eigenfunction. We have now 
obtained explicit expressions for or^_^ when k < 125 and oc^ m for m < 23 with arbitrary 
k. In order to do this the underlying symmetries of thecoefficients have to be exploited 
efficiently Some of this work was already done in Ref. [H|, and one can find explicit 
values of terms up to k < 49 and m < 10 in the addendum. 

We proceed by noting that the innermost coefficients, {o4-fc m }; yield a general 
solution to the present difference equation. We define 

^ J ( ;ir , n (»-")■ (it) 

and assume that all coefficients up to level a^\ k _-^ m are known. The requirement that 
the coefficient for each power of 5 cancels separately and some algebra yields explicit 
expressions for the coefficients on the next level. For even values of k we find 



a 



k'-i 

n-2l,l+m-k'l 11 ~t~ ^n,2Z+l,fc-2i«„_2i-l,i+m-A:'+l 



n—k,m>l 




l=k'— m— 1 



l=k'—m 

Correpondingly for odd values of k this yields 




a 



fe'-i 

( n ) I \ (n) / 

2_ ^n,2l+l,k-2ia n _2l- lt l +m ^ k '/ n 



n—k,m>l 



l=k'—m 
k< 

» 



,2l,k+l—2lCt n —2l } l+m—k' / (C n , fcj i/n - C n>fc)0 ). (19) 



l=k'—m 



Here aj^ 1 ,] = 1 and impossible coefficients are taken to be zero. The recursive equations 
do yield the general solution, but the complexity of equations grows at an exponential 
rate. Thus, it has been necessary to use even higher order symmetries to reduce the 
number of equations and reach the present order k = 125. All calculations have been 
performed using the symbolic mathematical software Mathematica. 

By restricting the allowed values of r to the set {fc^}^ we transform the problem 
into an eigenvalue problem for an infinite tridiagonal matrix. By denoting u k := u^\k5), 
we see that the vector {u k } h % 1 is an eigenvector of the matrix 

H kk = 5/k, H ktk+1 = H k+ltk = l/2, A; = 1,2,... (20) 



corresponding to the eigenvalue fi n = y 1 + (5/n) 2 . In addition, the exponential part of 
the solution simplifies to 

exp(-M arsinh(5/n)/5) = (y/l + (5/n) 2 - 5/n) k . (21) 

Results from numerical diagonalisation agree with our results within numerical precision, 
as long as convergence can be reached. For 5 real and positive, normalised eigenvectors 
form an ortho normal basis of £ 2 . 
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The polynomial part of the eigenvectors define a discretised version of the 
corresponding associated Laguerre polynomials L\. It is not clear whether a three-term 
recursion relation exists for these polynomials. 

4. Solutions in spectral representation 

The polynomial character of the exponential part shown in Eq. (121 J) also contains both 
the eigenvalue fi and the discretisation parameter, which indicates that we should also 
examine the problem with respect to the spectral variable. It turns out that the problem 
at hand corresponds to a special case of the Pollaczek polynomials jHj- They satisfy the 
three-term recursion relation 

(j + a, b) = 2[(j + A + a)x + b]Pf(x\ a, b) - (j + 2A - tfP^x; a, b), (22) 

where j > 0, and initial conditions 

P*(x;a,b) = 1, P£(x; a, b) = 2(A + a)x + b. (23) 
From Eq. (j2*U|) we obtain the recursion relation 

(j + l)u j+1 = 2[(j + l)x - 8] Uj - (j + j(= k - 1) > 0, (24) 

and identify the parameters A = 1, a = and 6 = —5. The discrete spectrum agrees 
with the calculations performed above, i.e. 



x m = ^/l + 5 2 /(m + l) 2 , m = 0,l,.... (25) 

We have not yet studied the absolutely continuous spectrum in the range [—1, 1]. 
The explicit formula for Pollaczek polynomials reads 

PnVsS; » = ± (-A + i*WMA + «W)„- V (a -,,) (26) 

where x := cos6», $(6») = (acos^ + 6)/sin6» and (A) k = A(A + 1) • • - (A + k - 1). The 
orthogonality of Pollaczek polynomial is defined with respect to the interval [—1,1], 
where cos 6 and sin 6 are easily defined. Here, it would be tempting to use 

sin# = ±i5/(m + 1), (27) 



which yields simple terms to be inserted in Eq. (J26J). Nevertheless, the correct way to 
do this as well as the corresponding interpretation are not obvious to us. 

Below, we obtain a simpler way to express the polynomials Pj(x) := Pj{x\ 0, —5) 
for a;(> 1) within the discrete spectrum (|23jl . The results of the previous section show 
that we can write 

Pj(x m ) = (x m - 5/(m + l)y- m Qf(x m ), (28) 

where Q™ is a polynomial of degree m. By extracting the polynomials Qj, we can 
reconstruct the corresponding Pollaczek polynomials. The next step is to evaluate 
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polynomial relations with respect to the index j, i.e. express the coefficients of the 
polynomials in as functions of j. Thus, we write 



C? + i)E 



Q j (x m ) = (j + l)J2 




(29) 



(30) 



Both factors (3j t i and ■jmj appear to be quite complicated at first. Some general features 
can be gleaned out, but the breakthrough is achieved in three steps. First is the 
observation that 7^ = fyj and the next amounts to the symmetry fijj, = fiij. Finally, 
we find 



ft, 



min(j,m) 




(31) 



This is the solution we have been looking for, explicitly: 

JM*m) = 0' + l)I> 
1=0 




= (./ + 1) x m - 

\ m + 1/ ^ \m + 1 

The first expression is more convenient for j < m, while the second is more compact for 
j > m. Note that expressions are manifestly identical for j = m. The unnormalised, 
general solution to the matrix eigenvalue problem (J2*U|) now simplifies to 

Wfc = Pfc_i(x„_i), fi n = X n -i. (33) 

For real values of 5 the solutions satisfy the orthogonality relation 



fc=i 



and normalisation requires that 



E 

fe=i 



u 



n\2 



1. 



(34) 



(35) 



For small values of j or m we can be sure that Pollaczek polynomial do satisfy the 
recursion relation. The solutions for spectral and coordinate representations are identical 
because they are solutions to the same difference equation. Thus, in the limit 5 — > 0, 
the Pollaczek polynomials tend to Laguerre polynomials with corresponding exponential 
parts, even if the variables x and r do not coincide. The extraction process guarantees 
that each term has been uniquely and correctly identified with corresponding powers of 
x m and —5/(m+ 1). Next, we identified several sets of equations that these coefficients 
satisfy. This allowed us to construct further coefficients in the series without redoing 
the extraction process. As the final step, the explicit expressions were conjectured and 
verified against known results. The general solution to Eq. (j2J) can be considered as a 
partially proven conjecture.^ 

I The required intermediate steps are available at request from Matias.Aunola@pvtt.mil.fi 
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5. Discussion 

We have examined a simple discretisation of the radial Schrddinger equation and shown 
that it is exactly solvable. We were able to obtain an explicit solution both in terms of 
the radial coordinate as well as the spectral variable, i.e. the eigenvalue. In future, one 
can simply use the existing solution, e.g. in the form of a piece of computer code based 
on Eq. (|32|). Simultaneously, we derived simple expressions to Pollaczek polynomials 
Pj(x] 0, —5) for the discrete mass points x m . Initial steps of the present approach are 
due to earlier research on the discretised ID harmonic oscillator, where we obtained 
asymptotical representations of Mathieu functions ^Hj. Much work still remains and 
alternative approaches should be applied to these problems. 

As a final note, we state that the present discretisation of the Schrddinger equation 
can also be used when visualising hydrogen radial wave functions. The possibility of 
comparing a numerical algorithm against exact results is not too common, especially 
if exact results are for the algorithm itself. Of course, it is not possible to carry the 
recursion either to infinite order or with infinite precision. 
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